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ABSTRACT 

During the epoch of reionization a large number of photons were produced with frequencies below 
the hydrogen Lyman limit. After redshifting into the closest resonance, these photons underwent 
multiple scatterings with atoms. We examine the effect of these scatterings on the temperature of 
the neutral intergalactic medium (IGM). Continuum photons, emitted between the Lya and Ly7 
frequencies, heat the gas after being redshifted into the H Lya or D Ly/3 resonance. By contrast, 
photons emitted between the Ly7 and Ly-limit frequencies, produce effective cooling of the gas. Prior 
to reionization, the equilibrium temperature of ~ 100 K for hydrogen and helium atoms is set by these 
two competing processes. At the same time, Ly/3 resonance photons thermally decouple deuterium 
from other species, raising its temperature as high as 10^ K. Our results have important consequences 
for the cosmic 21-cm background and the entropy floor of the early IGM which can affect star formation 
and reionization. 
Subject headings: intergalactic medium - diffuse radiation - radiation mechanisms: general 



1. INTRODUCTION 

During the epoch of reionization, UV photons between 
the Lya and Ly-limit frequencies are produced in abun- 
dance. Unlike the ionizing photons above the Ly-limit, 
whose mean free path in the neutral IGM is quite short, 
photons between Lya and the Ly-limit scatter only af- 
ter redshifting into one of the atomic resonances, which 
allows them to move freely through cosmological dis- 
tances. Eventually scattering these photons exchange 
energy with atoms, thereby providing a mechanism for 
raising or lowering the gas temperature far from any ra- 
diation source. 

Several authors previously considered this process. 
Madau, Meiksin & Rees (1997) estimated that photons 
lose about one per cent of their energy to the gas as they 
redshift through Lya resonance. They found accordingly 
that, if Lya scattering controls the HI hyperfine level 
population, as required to produce a 21-cm background 
distinct from the Cosmic Microwave Background (CMB), 
recoil heating sufficed to raise Tjom above Tcmb and 
to move the 21-cm background into emission. Later, 
however, Chen & Miralda-Escude (2004) showed that 
the inclusion of atomic thermal motion (neglected by 
Madau et al.) in the analysis reduces the heating rate by 
at least three orders of magnitude. Furthermore, they 
showed, while photons emitted between Lya and Ly/3 
frequencies (called by them the "continuum photons") 
heat the gas when redshifted into Lya resonance, those 
injected directly into Lya by the cascade from higher 
resonances (called the "injected photons"), cool the gas. 
With increasing gas temperature, the efficiency of "con- 
tinuum photons" as heaters declines, while that of "in- 
jected photons" as coolers rises. Therefore, if, as these 
authors claimed, the numbers of these two types of pho- 
tons were similar, then cooling would prevail at tempera- 
tures above 10 K. If so, the 21-cm background produced 
by Lya pumping would have been in absorption, unless 
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some other mechanism such as photoelectric heating by 
a strong X-ray background, can raise the IGM tempera- 
ture above Tcmb without fully ionizing it. 

However, two important physical mechanisms were 
missing in previous calculations. The first is the forbid- 
den transition from the 2s to Is level in hydrogen. The 
cascade which follows absorption of most photons which 
redshift into high atomic resonances (Ly/3, Ly7, LyJ etc.) 
proceeds preferentially not via the 2p level (whose decay 
to Is produces a Lya "injected photon") but via the 2s 
level, instead. This significantly reduces the number of 
injected photons, raising the equilibrium temperature to 
about 100 K. The second is the impact of deuterium. 
Despite its very low abundance, deuterium makes an im- 
portant contribution to the gas heating by its interaction 
with Ly/3 photons. At high radiation intensities, the D 
Ly/3 resonance becomes more important for heating the 
gas than H Lya. In addition, scatterings by D atoms re- 
duce the number of injected photons for hydrogen, thus 
decreasing the cooling rate. 

In §2 and §3 we calculate the effect of resonant scat- 
tering on H and D atoms, respectively. In §4 and §5 
we discuss the impact of this scattering on the thermal 
evolution, 21-cm signal and reionization history of the 
IGM. 

2. HYDROGEN RESONANCES 

When photons redshift through the Lya resonance, 
multiple scatterings affect their spectrum. The in- 
tensity J(i^) varies with frequency in the neighbor- 
hood of the resonance at Vg according to the following 
(|Chuzhov fc Shapiro 2001 

J{x) = J{Q)e-'^-'^\ (1) 

for injected photons and x > 0, otherwise 



J(x) = 27rJo7a ^e 3" 
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where a = A2i{2kTii/mc'^)~^^^ /iTTiyg, T^ is the hy- 
drogen kinetic temperature, A21 is the Einstein spon- 
taneous emission coefficient of the Lya transition, x — 



{v / Ua — V) / {2kTYi/rac^Y''^ is the frequency distance from 
line center divided by the thermal width of the reso- 
nance, 7"Hl + 0.4/Ts)~^ = TGP > 1 is the Gunn- 
Peterson optical depth, Jq is the UV intensity far away 
from the resonance. The recoil parameter, -q, equals 
[/ir/„/(2fcTHmc2)i/2][(i + o.4/r,)/(l + 0.4/TH)], where T, 
is the spin temperature of the 21 cm hyperfine transi- 
tion, which equals Tjj (Tcmb) at high (low) radiation 
intensities. The intensity at a; = 0, J(0), is given by 



j(0) _ <(J^(C)-J-^(C)) 



Jo 



V3 
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where C — {16rj'^a/9TTj)^''^, 1F2 is a hypergeometric func- 
tion and Ji and J_ 1 are the Bessel functions of the first 
3 3 

kind. 

Since photons to the red (blue) of the resonance scat- 
ter preferentially off atoms moving towards (away) from 
them, the average energy an atom gains from a scattering 
photon depends on its frequency: 



AEi.) = i^ (1 



h <j){x) 



(4) 



where <ii(x) is the norma lized scattering cross-section 
(jChuzhov fc Shapiro 20051 ) . The total energy gain of the 
gas from each photon as it passes through the resonance 
is determined by the radiation intensity, the probability 
of scattering and the average gain at each frequency: 
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(5) 



Figure[l]shows the heating and the cooling rates obtained 
by numerical integration of equation ^.^ For Th > 100 
K and the range of optical depths corresponding to the 
mean density IGM at 10 < 2; < 30, AE'tot can be well 
approximated by the asymptotes 



AEjk - 0.37Th ^/^(l + z) K 



(6) 



for continuum photons and 

AEi/k^-Q.3T^^^{l + z)^^^K (7) 

for injected photons. The total heating/cooling rate is 



H^ = Mae, + ^AEi) 



(8) 



where Na is the number of the photons that pass through 
the Lya resonance per H atom per unit time. 

The continuum photons (i.e., those that are red- 
shifted into the Lya resonance) can be produced in 
two ways. Most originate as photons emitted be- 
tween Lya and Ly/3, which eventually redshift into the 
Lya resonance. The rest come from higher-frequency 
photons that are absorbed by D atoms and cascade 
into the D Lya resonance, from which they then red- 
shift into the H Lya resonance. The injected pho- 
tons (i.e., those injected directly into the Lya reso- 
nance) come from the photons emitted between Ly7 and 
the Ly-limit. Most of the latter photons, when red- 
shifted into the closest resonance and absorbed by H 

^ Results are virtually insensitive to the choice of Lorentz vs 
Voigt profile for <j>{x). 



atoms, produce a cascade to the 2s level, from which 
an electron decays to t he ground level by emitting 
two photons below Lya (las also noted bv Hirata 2006t 
iPritchard fc Furlanetto 2006f ). Another fraction are ab- 
sorbed by D atoms and either destroyed or converted to 
D Lya photons, as previously explained. The rest, after 
redshifting into the closest H resonance and producing 
a cascade to the 2p level, become the injected photons. 
Therefore, the ratio of injected and continuum photons. 
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where Vn is the frequency of np- 



t'lim(l 



'), where i/ui, 



• Is transition (i.e., 
13.6 h-^ eV is the 



Ly-limit frequency), P„ is the fraction of all cascades 
from level np that go through 2p (0, 0.26 for n — "i and 
4, respectively, and roughly 1/3 for n > 4), Js{v) is the 
spectral profile of the radiation source, and /D,n is the ra- 
tio of the cascades from level np occurring in deuterium 
and hydrogen. It can be shown that, since the deuterium 
optical depth, rD,n, in the np —t Is transition, was low 
during reionization, /d.ii ~ 1 — e-'rD.n(i-/rcc,n)^ where 
/roc,n is the probability for the electron in level np to 
decay directly to Is. For a hot source with surface tem- 
perature Tbb ^ 10^ K, Jsiy) (X v^ between Lya and Ly- 
limit, and Ji/Jc ~ 0.17 (the exact value depends on the 
optical depth of deuterium which is redshift dependent). 
Colder sources would produce a lower Ji/Jc ratio, but 
unless the spectrum color temperature is below ^ 5 • 10** 

K, Ji/Je > 0.1. 

Combining equations ([6]) and (l7|), we find that, for the 
mean IGM, the equilibrium temperature at which cool- 
ing by injected photons balances heating by continuum 
photons is 

1/2 / 7. / T \ -3/2 



Teq « 130 K 



10 



Ji/Jc 

0.15 



(10) 



This greatly exceeds the mean temperature prior to 
reionization, if only adiabatic cooling occurs once the 
IGM decouples from the CMB. 

3. DEUTERIUM RESONANCES 



Unlike hydrogen, the 
nance for deuterium is 



most important res- 
not Lya, but Ly/3 
As photons redshift 



(Chuz hov fc Sha piro 2005). 
through the D Ly/3 resonance, a significant fraction of 
them are destroyed via absorption and cascade. Prior to 
reionization (z > 10), the photon destruction probability 
is above 0.2. Therefore the blue wing of the resonance 
is significantly higher than the red wing, so that the 
radiation color temperature around the resonance is 
negative. Because of their negative color temperature, 
the D Ly/3 photons are much more efficient than the 
Lya photons at heating the gas. 

If the velocity distribution of deuterium atoms is 
Maxwellian, then the radiation intensity around the Ly/3 
resonance varies according to 

J(x) - 7d J'(x) = 



^ 



0.5Acc.3 / ETic{Max[\x'\,\x\])J{x')dx', (11) 



where "■ 



7r.- = Tn.., ^ 3 • 103(1 + z)3/2(„D/nH) 

Chuzhoy &: Shapiro 20051 ). Solving equations ([5]) and 



(|lip numerically, we find that for neutral IGM between 
z = 10 and 20 the average energy each Ly f3 photon trans- 
fers to the gas is 



AEtot/k « 0.012Ti 



1/2 



1 + z 
15 



3/2 



K, 



(12) 



where Td is the deuterium kinetic temperature. 

Physical insight into the above result can be gained 
from the following argument. Photons are scattered by 
atoms in whose rest frame they are close to resonance 
(i.e., iy{l — v/c) ~ ^'Ly/3, where v is the velocity of an 
atom). Thus photons in the red (blue) wing of the res- 
onance are scattered preferentially by atoms moving to- 
wards (away) from them. When the spectrum around 
the resonance is flat, then to first approximation, the 
energy gain of atoms from photons in the blue wing is 
compensated by energy loss to photons on the red wing. 
However, when the spectrum is tilted, as in our case, 
each extra photon on the blue wing adds to the gas 

1 /2 

AEtot/k -- {hi^Lyp/k)vthermai/c w 0.05Tj^' K, Compa- 
rable to what we got in equation p^ . We note that, 
unlike H Lya photons, D Ly/3 photons become more ef- 
ficient heaters when the temperature rises. 

Since radiation affects H and D atoms differently, we 
now have to make separate estimates of their respective 
temperatures ^. From equation (|12p we find that the 
heating rate per D atom is 



10-2« erg • s-i 



2-10-5 




(13) 



where N/j is the number of photons passing through the 
Ly/3 resonance per H atom per second. Elastic collisions 
with H atoms and to a smaller extent, with He atoms, 
make D atoms lose energy at a rate 
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where cthd is the momentum transfer cross-section for 
collisions between H and D atoms "'. When the radiation 
intensity is high, so that Td ^ Th, Td is set by the 

^ In principle the same is also true for other species (He, p, e), 
but in practice only the temperature of deuterium is significantly 
different. 

* Some element of uncertainty is present in the above equation 
due to the absence of experimental measurements of cthd- The 
collisional cross sections of atoms belonging to different species 
can be roughly determined by the size of t he atom, which for D 
and H atom s gives (Thd = 0.35 • IQ— '^^cm^ (IClementi et al. 19631 : 
ISlater 1964t ). However, since D and H atoms energy levels are 
nearly degenerate it is possible that cthd is significantly larger. 




Fig. 1. — The heating/cooling rate due to continuum (solid 
and dash-dotted lines for Ts = Th and Ts = Tqmb respectively) 
and injected (dashed line) H Ly« photons, and D Ly/3 photons 
(dotted line) for fiat spectrum source (Ji/Jc ~ 0.15, J^/Jc ~ 0.35) 
and 2 = 12. For Lyo and Ly/3 photons the results are plotted vs 
temperature of H and D atoms, respectively. 



balance between heating and cooling rates, Xd = H^ 
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4. THERMAL EVOLUTION OF IGM 



When resonance photons are the only important heat- 
ing source, the temperature of H and D atoms evolves 
according to 



dt "3fc^ " 
dT^_ 2 
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(16) 
(17) 



where the term 4T/3i accounts for adiabatic cooling due 
to Hubble expansion. The hydrogen temperature at any 
redshift, rH(^), does not depend very strongly on the his- 
tory of photon production (see Fig. [2]). Thus Th(z) may 
be expressed as a function of the total number of photons 
between Lya and the Ly-limit per baryon emitted by red- 
shift z, N^^Yi (Fig. El- If Pop HI stars are the main radia- 
tion sources, then, by th e end of reionization N-^^yi ^ 200 
(jCiardi fc Madau 20031 ). This implies that Th rises to 
~ 100 K, while Td is raised to - 10'*(cthd/10"^5 cm^y^ 
K. However, both much higher and much lower maximum 
values for N^^-a are not completely excluded. 

5. IMPLICATIONS 

We have shown that, during the reionization epoch, 
resonance photons can heat the neutral IGM to about 
lOO K, but not much higher. When gas collapses out of 
the IGM, the entropy cannot decrease (without radiative 
cooling), so its overdensity inside virialized halos cannot 
exceed (Tvir/To)^''^, where Tq (Tvir) are the temperatures 
before (after) collapse, respectively. Consequently, for 
minihalos with Tvir < 4000 K that form out of gas pre- 
heated to 100 K, the post-collapse overdensity is less than 




Fig. 2. — Thermal evolution of H (solid lines) and D (dashed 
lines) atoms for cthd = lO-^^ cm^ and Nj^h(z = H) = 200. The 
comoving radiation intensity is assumed to grow as e~^' ^^, where 
Az = 1, 2, 3 correspond to lower, middle and upper lines. 
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Fig. 3. — Hydrogen kinetic temperature at 2 = 15 (solid and 
dotted lines) and 2 = 11 (dashed and dashed-dotted lines) versus 
^7,ij(-z), for Cud = 10~^^ cm-^. The dotted and dashed-dotted 
lines are the results when the heating of deuterium is neglected. 

~ 200. Minihalos with such lowered central densities 
form H2 molecules more slo wly and are more vu lnera- 
ble to H2 photodissociation (|0h fc Haiman 20^31 ). The 
photodissociation is caused by UV photons in the range 
from 11.6 to 13.6 eV, which overlaps the range from 10.2 
to 12.7 eV (from Lya to Ly7) responsible for preheating. 
This may suppress H2 cooling and star formation in these 



minihalos. Minihalos would also be more easily photoe- 
vaporated, thereby reducing their consumption of ioniz- 
ing photons during reionization (Shapiro, Iliev & Raga 
2004). 

This preheating also has a strong impact on the red- 
shifted 21-cm signal from the beginning of reionization. 
UV resonance photons may raise Th above Tcmb, thus 
moving the gas from absorption to emission (see Fig. |4]). 
In fact, even relatively low intensities (< 10 photons per 
baryon) can drastically affect the signal. 

If the first radiation sources deposit a substantial frac- 
tion of their energy in X-rays, th e average temperature of 
the IGM may climb even higher ( iVenkatesan et al. 20011 : 
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Fig. 4. — The evolution of H 21-cm differential brightness tem- 
perature for N^^h{z = 11) = 10 with (solid) and without (dotted) 
taking the gas heating into account, and for N-y^ni^ = H) = 200 
with (dashed) and without (dashed-dotted) gas heating. The 
comoving radiation intensity is assumed to grow as e~^'^ and 
o-HD = 10" 



-15 



lOh fc Haiman 2003f ). However, unless the X-ray spec- 
trum is very hard, most of the energy is deposited within 
~ 1 comoving Mpc of the source. By contrast, the UV 
photons considered here, which travel cosmological dis- 
tances before scattering, produce a very homogeneous 
heating mechanism. Therefore, in regions that do not 
have a radiation source nearby, heating by resonance 
photons can be more important even if, on average, X- 
ray heating is stronger. 
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